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Abstract 

Euclidean invariant Klein-Gordon, Dirac and massive Cher n- Simons 
field theories are constructed in terms of a random walk with a spin factor 
on a three dimensional lattice. We exactly calculate the free energy and 
the correlation functions which allow us to obtain the critical diffusion 
constant and associated critical exponents. It is pointed out that these 
critical exponents do not satisfy the hyper-scaling relation but the scaling 
inequalities. We take the continuum limit of this theory on the basis of 
these analyses. We check the universality of obtained results on other 
lattice structure such as triclinic lattice and body centered lattice. 
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1. Introduction 



The random walk method has been studied extensively as a useful tool in sta- 
tistical mechanics as well as in constructive field theory. Many critical phenomena 
have been revealed in this language. For instance, in the two-dimensional Ising 
model, the free energy and the correlation functions are evaluated by summing over 
all diagrams in the high temperature expansion in terms of random walks with a 
spin factor [|IJ. This study also shows that Majorana spinor field can describe crit- 
ical phenomena. As for ip 4 theory, the studies of Brownian motion on an arbitrary 
dimensional lattice enables us to prove important correlation inequalities in terms 
of a random walk representation 0. The construction in three-dimensions and the 
proof of triviality in higher than four dimensions can be done using this method. 
Although there have been many rigorous studies on constructing a spinor field by 
random walk method on a two dimensional lattice and also on Brownian motion on 
an arbitrary dimensional lattice for scalar field, a random walk construction of a 
spinor field on a higher dimensional lattice has so far never been obtained. 



The purpose of this paper is to construct a theory of free spinor fields by random 
walks on a three dimensional lattice. We employ a regularization such that all quan- 
tities can be evaluated as finite ones at any stage before taking the limit of lattice 
spacing to be zero. So far there have been some studies of a continuum random walk 
method for spinor fields ||-||. Consider the following two point function of spinor 
fields given in terms of a path integral representation pfl: 



< x, 



< e, 



e,; > <g> X; > = 



m 



Dx 



exp 



Diff 



-m 



Jo [ V¥ l 

(1.1) 

where J M is spin matrix with the magnitude J and exp[zJ$] is the spin factor. $[e] 
is the oriented area determined by path e(t) between = e(0) and = e(l) on a 
two dimensional unit sphere ||. The boundary condition of the path integral Ql.lD 
requires initial and final data 



X; 



The state vectors < e^| and |e^ > are SU(2) coherent states for the spin matrix 
J M In this paper, we give a mathematically well-defined representation of the 
correlation function ( |1 . 1|) with a lattice regularization. 
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First we work on a three dimensional regular cubic lattice. We define a quantum 
diffusion process in terms of an amplitude with a spin factor whose discrete time 
evolution obeys a certain recursion relation. This random walk model is a simple 
generalization of the diffusion process for the two dimensional Ising model to the 
three dimensional model. The free energy and the two point function are calculated 
by solving the recursion relation. We show that there are only J = 0, 1/2, 1,3/2 
and 2 cases on the regular cubic lattice. The critical diffusion constant for each J is 
found at the singular point of the free energy and its associated critical exponents 
are obtained. This analysis shows how we should scale the variables to obtain the 
two point correlation function of the continuum field theory in our random walk 
model. 

We obtain Euclidean invariant field theories with spin 0, 1/2 and 1 for J = 0, 1/2 
and 1, respectively. However, theories with spins higher than 1 are not obtained from 
this construction. This is an unexpected result which cannot be understood in a 
naive continuum random walk method for spinor fields. 

The classification of universality in J = 0, 1/2, 1 models is nontrivial. The order 
of phase transition of the model with J = and 1/2 is shown to be second and third, 
respectively. In the J = case, the model belongs to the same universality class as 
Brownian motion whose typical path has Hausdorff dimension 2. The Klein-Gordon 
field is constructed at this critical point. In J = 1/2 and 1 typical path has 

Hausdorff dimension 1. The free Dirac field is constructed at the critical point of 
J = 1/2 model. It is worth noting that there is no fermion doubler in this model. In 
the J = 1 case, the critical behavior looks strange due to the existence of Euclidean 
non-invariant excitations with long range correlations in addition to an invariant one. 
They contribute to the most singular terms of the free energy, which suggests that 
the order of phase transition is first. Nevertheless, we can take a suitable continuum 
limit so that they do not appear in the correlation function. The massive Chern- 
Simons field theory is constructed at this critical point of J = 1. These non-invariant 
massless excitations affect critical exponents which do not satisfy the hyper-scaling 
equality dv = 2 — a but satisfies the inequality dv > 2 — a. 

In the J = 3/2 rotational invariant correlation function is never obtained. 

In the J = 2 case, there are two critical points for real k, one of them gives an 
Euclidean invariant scalar field theory but the other does not. 

We also investigate the lattice dependence of our model. The universality of 
our construction should be checked on other lattice structures. We examine it on a 



3 



triclinic lattice with bonds having identical length and on a body centered lattice. 
The results show that the spin and the spin 1/2 field theories are also obtained on 
these lattices while the spin 1 field is not obtained in general. 



This paper is organized as follows. In section 2, we give the definition of the 
random walk model on a three dimensional cubic lattice. We study the critical 
behavior of this random walk model in section 3. In section 4, the Klein-Gordon, 
Dirac and massive Chern-Simons fields are constructed by taking a scaling limit. 
Various exponents are shown for convenience. In section 5, we attempt to check 
the lattice universality of our results obtained in the previous sections, especially 
on triclinic and body centered lattices. Some remarks and comments are given in 
section 6. 



In this section, we consider a random walk model on a three dimensional cubic 
lattice Z 3 , which is a discrete version of fll.l]) 



where J2u> denotes the sum over all closed non-backtracking walks on the cubic 
lattice, k is a diffusion constant and J is a half-integer valued parameter which will 
be interpreted later as spin of the model. The quantity S(u) is the oriented area 
encircled by unit tangent vectors of a walk u on S 2 . The weight for a path in this 
model is not positive-real- valued but complex-valued. It causes cancellation between 
walks, and the mean free distance of random walks becomes longer than that for 
Brownian motion. Such smoothing affects the critical behavior of our model. 

Let us consider cubic lattice generated by three orthonormal vectors e!,e2,e 3 
which satisfy e a ■ ep = 5 a p (a, (3 = ±1,±2,±3). Here and after we use notation 
e_ Q = — e a . This means the lattice unit is taken to be unity. A site x is represented 
by x a G Z as x = J2a=i x a e a . In order to sum up all possible walks, it is convenient 
to introduce an amplitude Af (x|«, J), which is the sum of contributions from all 
walks which start from the origin with the direction e a and arrive at a site x from 
the direction after /-steps Fig.|l|. 



2. 3D Random Walks with Spin Factor 




(2.1) 
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Figure 1: Figure of a walk start from the origin with the direction e Q and come into 
a site x from the direction eg and the cooedinate system. 

As we will see below, (|2.1|) can be represented on a lattice as a memory- 1 walk. 
The weight for a walk can be divided into two parts. One is fi>l/H which depends 
only on the length of the walk uj. The other is the spin factor e lJS ( u ) which gives 
non-zero phase when a walk curves. To define such a phase on the lattice, we 
need successive two steps of a discretized walk. The amplitude obeys the following 
recursion relation 

Ag^xlK, J) = k e iJ5(e3 ' e " >e ^ 7 (x - e/9 |«, J), (2.2) 

where —a denote the direction of — e a and S(e a , ep, e 7 ) is the area of a spherical 
triangle with vertices e a , ep and e 7 on S 2 , basically. The precise definition of the 
phase factor however, will be given later. The sum is taken over all nearest neighbour 
bonds except the direction — (3 in order to forbid backtracking processes. The initial 
condition is given as 

Af(x.\ K ,J) = <F% l0 <W*sO. (2-3) 

Now we define the free energy density (free energy per unit volume) in terms of 
the amplitude as 

/(*, J) = -(-1) 2J hm E Ur(0\K, J). (2.4) 
L ^°° i=i Q t±i 1 

The factor l/l in this definition removes I degenerate walks which represent the 
same configuration in the sum, as in the case of two dimensional Ising model. The 
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factor (— 1) 2J attaches a minus-sign to every fermionic loop. This is important to 
construct both bosonic and fermionic theories in the same framework. 



The free energy density eq^^) gives an explicit definition of eq.( [2.1|) . The 
Fourier transformation of the amplitude A" l3 (p\n, J) = ^^A^^^k, J) e * p x enables 
us to solve the recursion relation 



ir+i(pi«> J ) 



k jr Ar(p\K,J)Q^(p,J). 

7=±1 



;2.5) 



The matrix Q a/3 (p, J) is a propagation matrix. The elements of the matrix Q af3 (p, J) 
are given as follows: 

Q(p,J) = Q(0,J).P(p) 



Q(0,J) 



-ILi J 

e 2 

e -2niJ 



Hi J 

1 



-i J 
1 

— ILi J 

e 2 

1 



3 27ri J 
g7ri J 

1 
1 

,—iriJ 





V 







-vp\ 



1 
1 

— ILi J 

e a 



1 





-IP2 



--H 



e niJ 
ILi T 

1 
1 



1\ 
1 


1 
1 



(2.6) 



By ( |2.6| ), we give the definition of the spin factor. Several comments must be noted. 
The representation of the matrix Q(0, J) is not the unique choice to obtain the 
same free energy density as ( |2.4| ). The zero components in the propagation matrix 
show the non-backtracking condition. The parameter J must be chosen to be a 
half-integer 

J G Z/2 

in order that the total phase factor for any closed path produces oriented area on 
S 2 . Note that the matrix Q al3 (p, J) is hermitian at p = 0. 

Using eq.( [2.5|) and the condition (|2.3|) the free energy density is reduced to the 
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Figure 2: An assignment of phase-factor. 



following form 

fM = -(-i)" & £ r ^ s tr|i[«g( Pl j)]« 

= (-l) 2J £^trlog[l- K g(p,J)] 

= (-l) 2J £^logdet[l-«;Q(p,J)]. (2.7) 

Here the second equality of the above equation is only formal one, because the 
convergence of the infinite series of the r.h.s. of the first line in eq. fl2.7Q is not 
ensured. In appendix, we show that there exists a critical diffusion constant k c such 
that for any |ac| < |/c c | the infinite summation converges uniformly with respect to 
p. There it is shown k c is the inverse of the eigenvalue of the matrix Q(0, J) whose 
absolute value is not less than that of other eigenvalues. The free energy is singular 
at k — k c . As we will see, a continuum field theory is constructed by taking the 
limit k — ► k c . 

Therefore we have to calculate the eigenvalues of the matrix Q(0, J), that is, to 
solve the equation 

det[x - Q(0, J)] = 0. (2.8) 
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Note that this equation has the following symmetries 



J -> J + 4, 

J -> -J. (2.9) 



The former is obvious from the eq. (|2.6| ). The latter is due to the fact that the 
transposed matrix of Q(0, J) is equal to Q(0, —J)- Moreover, we can show that the 
free energy possesses these symmetries as well. Thus, there are only five independent 
models corresponding to 

j={o, i, 1, | 2}. (2.10) 



By solving the algebraic equation (|2.8j) of six-degree in x, we can obtain all 
eigenvalues of Q(0, J). They are listed in table 1. The degeneracy of the eigenvalue 
with the largest absolute value for each J is 1, 2, 3 and 4 for J = 0, 1/2, 1 and 3/2 
respectively. 

We represent the explicit form of the free energy density for each values of spin 

J. 

f(K,J) = (-1) 2J f -^-x/(p|«,J), 



J-* (27r) 3 

3 

/(p|K,J=0) = log(/« 2 -l)(K 2 + l)(5/s: 2 -2fi;^cosp i + l) 

i=l 

3 



/(p|«,J=l/2) = log 



(57 - 40V2> 6 + 3k 2 + 1 + 2(-5 + 4V2)k 5 ^ 

COS Pi 



1=1 



<Spm Factor 




J=0 


5 


1 




1 1 


-1 


-1 


J=l/2 


1 + 2^2 


H 


-2V2 


1 - V2 l-y/2 


1 - y/2 


1 - V2 


J=l 


3 


3 




3 -1 


-1 


-1 


J=3/2 


1 + V2 


H 


-V2 




1-2V2 


1 - 2^2 


J=2 


3 


3 




1 1 


1 


-3 



Table 1: Table of eigenvalues of the matrix Q(0, J). The eigenvalues with the largest absolute 
value are written in boldface. 
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/( P kJ=i; 



+k 4 ^3 + 8 ( — 1 + y/2) (cos pi cos P2 + cos P2 cos p 3 + cos p 3 cos pi ) 

/ 3 3 \ 3 

+4k 3 i^cospi + 2(2 - \/2) J|cospi ] - 2k ^ cos Pi 
log [-27k 6 - 9k 4 + 3k 2 + 1 



i=l 



i=l 



/( P |k,J = 3/2) = log 



3 / 3 3 

18k 5 ^ cos + 4k 3 cos + 4 cos Pi 

j=l \i=l i=l j 

3 

(57 + A0V2)k 6 + 3k 2 + 1 - (10 + Sv^k 5 ^ cos ^ 

1=1 

+ K 4 (3 — 4(1 + V / 2)(cOSpi COSP2 + COSP2 COSP3 + COSP3 cos pi 
/ 3 3 \ 3 

+k 3 4 ^ cos pi + 8(2 + y/2) J[ cos pi ) - 2k Y cosp 

\ i=l i=l J i=l 



(2.11) 



/(p|«,J=2) = log 



-27k 6 + 3k 2 + 1 + 30k 5 

COS Pi 



1=1 



-k 4 (9 + 32(cospi cosp 2 + cos P2 cosp 3 + cosp 3 cospi)) 



+4k 3 Y, cos Pi + 8 Y[ cos Pi J — 2k cos Pi 



\i=l 



1=1 



i=l 



The two-point function is defined as 



iT a/3 (x|K, J) = lim yYQ-^O, J)) ai ■ Af(x\K, J). 



(2.12) 



1=0 



Here, the sum is taken over all walks which start the origin toward the direction e a 
and come into a site x from the direction ep. Using eq.fl2.5|) and the condition ( |2.3| ), 
the two-point function (|2.12|) is rewritten as 



^(x|k,J)= £^e ip - x ^(p|K,J), 



(2.13) 



where we use the following notation 



tf°"(p|«, J) = [[(1 - kQ(p, J))Q(0, J)] 



-1 



a/3 
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3. Critical Behavior 



In the previous section, we have defined a random walk model with complex- 
valued spin factor. In this section, we argue critical behavior of our model near k c 
< \k c \). As we have already discussed, there are only five independent models 
on the cubic lattice. They correspond to the case J = 0, 1/2, 1, 3/2, 2. We discuss 
the critical behavior for each J. Here we evaluate the singular part of the free energy 
(|2.7| ) with respect to the following variable 

s = . 

K c 

1) J = case 

In the case of J = 0, the phase factor is trivial as shown in the eq. ( |2.6| ). This model 
is well-known as Brownian motion. The free energy density at low energy is 



/(k,J = 0) 



p~° d 3 p 



(27T) 



log 



576 
3125 



(p 2 + 4s) + 



(3-1) 



The notation p ~ denotes the expansion around p = 0. The most singular part 
of the free energy density with respect to s is 



/(k, J = 0) ~ s^logs (k/1/5). 



(3.2) 



Here we see the phase transition is of second order, therefore the exponent a = | is 
obtained. This result agrees with the well-known result of Brownian motion. 

2) J = 1/2 case 
The low energy limit of the free energy density for J = 1/2 is 



/(«,j=i/2)= r°^\o< A 



36 



(2tt) 3 to L(l + 2 v / 2) 4 



p 2 + 9s 2 ) + 



The most singular part of the free energy density is 

/(«, J= 1/2) ~s 3 logs (k/ (1 + 2V/2)" 1 ). 



(3.3) 



(3.4) 



Therefore the phase transition is of third order, and the exponent a = — 1 is ob- 
tained. 
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3) J = 1 case 

The critical behavior of the J = 1 model is more complicated than the former two 
cases because there are several massless excitations. The singular part of the free 
energy is dominated by the integral region near p = (tt, 0, 0), (0, tt, 0), (0, 0, tt) as 
well as p = : 



/(«,■/=!) = / 



p~o ftp 



(27T) 



log 



+ t(PiW + P2 W + VzV\ 



+ 



+ 



p~M,0) 


d 2 p 


log 


■16 




(27T) 2 


.27 


P~(0,7T,0) 


d 2 p 


log 


■16 






.27 


P~(0,0,tt) 


d 2 p 


log 


16 




(2vr) 2 


.27 



i(pi 2 +p 3 2 + 4s) + 
i(pi 2 +p 2 2 + 4s) + 



(3.5) 



The last three excitations do not recover the three-dimensional rotational invariance 
even near the critical point. They are all two-dimensional scalar fields. The most 
singular part in s of the free energy density comes from the last three terms in ( fj.5| ) 



/(«,J = l)~slogs (k/1/3). 



(3.6) 



We can read from this equation that the phase transition is of first order *. 
4) J = 3/2 case 

In this case, there are three excitations near p = (71% 0,0), (0,7r,0), (0,0, ir) as well 
as in the J = 1 case. The degeneracy of the eigenvalue 1 + \[2 of the matrix 
Q(p, J = 3/2) for p = is 4, but that of the three other excitations is 2. The low 
energy scheme of the free energy density for J = 3/2 leads to 



/(«, J = 3/2) 



p~0 (ftp 



(27T) 



log 



32 



J (Pl 2 P2 2 + V2V 2 + P3 2 Pl 2 ) ) 



*One must exercise care to determine the true order of the phase transition for J = 1, because 
the excitations near p = (ir, 0, 0), (0, n, 0), (0, 0, ir) resemble those in two dimensional Brownian 
motion which is ill-defined. In practice, two dimensional Brownian motion has I.R. divergences. 



11 



+ J 


^•P~(7T,0,0) 


d 2 p 

(27T) 2 


log 


+ J 


f P ~(0,7r,0) 


d 2 p 
(2tt) 2 


log 


+ J 


fp~(0,0,i") 


d 2 p 

(27T) 2 


log 



32 



(1 + V2) 2 
32 

(1 + V2) 2 
32 

(1 + V2f 



(P2 2 + P3 2 + 4S 2 ) + 

{Pi 2 + Ps 2 + 4s 2 ) + 
(Pi 2 +p 2 2 + 4s 2 ) + 



(3.7) 

Note that the term p 2 p 2 +P2 2 pz 2 +Pz 2 Pi 2 breaks rotational invariance. The most 
singular part with respect to s is 

/(«, J = 3/2) ~ s 2 logs (k / (1 + v 7 ^) -1 ). (3.8) 

The phase transition is of second order, and the exponent a is 0. 
5) J = 2 case 

There are two critical points k c = 1/3 and k c = —1/3 in this case. In the case of 
k c = 1/3, the low energy limit of the free energy density is 



/(«~1/3,J = 2) = J 



p~° d 3 p 



(2vr)3 l0g 



— s(p + s) 

81 V V 2 ' 



+ ^(PiW + P2 V + P3 2 Pi 2 ) ) H 



(3.9) 



There exists a non-rotational invariant term. The leading behavior of the free energy 
density in terms of s is 

/(k, J = 2) ~ sf logs (k/1/3). (3.10) 

The phase transition is of second order, and the exponent a is 1/2. 

In the case of k — — 1/3, the low energy limit of the free energy density is 

"64 



/(«~-i/3,j=2)= r°^o g 



81 



(p 2 + 12s) + 



(3.11) 



The leading behavior of the free energy density in terms of s is 

/(k, J = 2) ~ sf log s (k \ -1/3). 
The phase transition is of second order, and the exponent a is 1/2. 



(3.12) 
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4. The Scaling Limit 



In this section, we derive the n-point correlation functions for continuum spinor 
field theories. 

We restore the lattice spacing parameter a by the following replacement : 

x/a 



x 

p -> pa 

s = (m • a) 



(4.1) 



Here the exponent £ must be chosen in such a way that the correlation functions have 
suitable continuum limit, the scaling limit. This can be achieved by taking a — > 0, 
s — > while keeping m finite. We investigate the scaling limit of the correlation 
functions. We have to evaluate only the two-point function since all the theories we 
construct are free field theories. We can see whether they have rotational invariance 
or not. 

To extract the massless mode from the two-point function K (xI3 {x\k 1 J), we choose 
a unitary matrix U(0, J) which diagonalizes the mass matrix 1 — kQ(0, J). For this 
purpose, we introduce 



/rsL(x|«, J) = u(o, j) • /^(x|«, j) ■ u(o, jy 



(4.2) 



We choose v = 1/2, 1, 1 for J = 0, 1/2, 1 respectively for the scaling limit. The 
dominant part of the two-point correlation functions for J = 0, 1/2, 1 becomes 



#2L(x|k,J = 0) = ^ 



4a /■*/« . 



/^l(x|«,J=l/2) = 



5 J-TT/a (2tt) 3 

7r/a 





0(1 5) 


p 2 + 4 m 2 


0(5 1) 


0(5 5) 



\ 



0(a 2 ) 



3a 2 



1 + 2-\/2 (2vr) 3 

+0(a 3 ) 
2a 2 / >7r / a d 3 p 





0(2 4) 


kr p + 3 m 1 (22) 


0(4 2) 


0(4 4) 



A diag ( X |«,j-ij - 3 y Wa(27r) , 
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0(3 3) 


iLp + 2 m I (3 3) 


0(3 3) 


0(3 3) 



+0(a 3 ). 



+ 



p indep. 
matrix 



(4.3) 



Here 0(ij) denotes an i x j zero matrix and La) is the i x % identity matrix. The 
matrices <T are the Pauli matrices and the 3x3 matrices L are the spin-1 matrices 
defined by [Li,Lj] = ie^L^ and L ■ L = 2. The p independent matrix appears in 
i£"di ag (p|/«, J = 1) is originated from the three other excitations which do not produce 
long-range correlation in the scaling limit. Actually, the p independent part gives 
a 5-function in coordinate space. We can subtract such a short range function from 
the two-point correlation function. Note that the rotational symmetry breaking 
terms in the two-point functions corresponding to those in eq. (|3.5| ) disappear in the 
scaling limit. 

In the J = 3/2 case, however, the situation is completely different from J = 1. 
We can not obtain the three-dimensional rotational invariant theory in the scaling 
limit. Actually, whatever the choice of v we cannot remove the breaking term 
Pi 2 P2 2 + P2 2 P3 2 + p 2 p 2 contrary to the J = 1 case. 

The correlation functions for long range modes are 



3(x|J = 0) 



def 



<7(x|J=l/2) 



def 



<7(x|J=l) 



def 



-a/3 



a->0 4a 

1,2 . n 2 



3 ipx 



d 3 p 

1 + 2^/2 



(x|k, J = 0) 






0(15) 


p 2 + 4 m 2 


0(5 1) 


0(5 5) 



lim 



K d t K (x|«,J=l/2) 



d 3 p 
-oo (2vr) 3 



3 ipx 



1 


0(2 4) 


icr-p + 3 m I (22) 


0(4 2) 


0(4 4) 



!S 2^ -i^W/c, j = i; 

s=m-a 
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- f 

J —I 



d 3 p 

J2^y 





0,3 3, 


iL-p+ 2 m I (3 3) 


0(3 3) 


0(3 3) 



(4.4) 



/ 



where we removed the 5-functions of x in the J = 1 case. We have constructed the 
two-point correlation functions for the Klein-Gordon, Dirac and the U(l) Chern- 
Simons fields from our random walk model with J = 0, 1/2 and 1 respectively. 
Therefore these free field theories are obtained rigorously. On the other hand, the 
spinor theories with spin higher than 1 cannot be obtained from this construction. 

Now we discuss critical exponents and relations among them. The susceptibility 
matrix can be defined by analogy with classical spin systems 



1 



a/3 



V(l-«Q(0,J))-Q(0,J)J diag (45) 

The susceptibility exponent (matrix) ■J is defined from the behavior of \ near the 
critical point in s 

X&~(«c-«)-T*~ a -7* (4.6) 

We see the matrix / ") /Q!/3 is the identity matrix of order 1, 2 and 3 for J = 0, 1/2 and 
1 respectively. 

The power behaviors in |x| at long-distance of the correlation functions are 

1 



lim g(x.\J = 0) 

x|^oo 

lim <7(x|J = l/2) 

|x|^oo 

lim g(pc\J = 1) 

|x|^oo 



|X| 

<r ■ e (x) 



|x 
L e 



2 

(x) 



12 ' 



(4.7) 



where e( x ) = (xi/|x|, x 2 /|x|, x 3 /|x|). Thus the exponent r] for each J is r\ = 0, 1, 1 
for J = 0, 1/2, 1 respectively. 

The exponential behavior at large-distance is also important. The correlation 
length is the inverse of the mass-matrix 



3£g (((i_«Q(o,J)).Q(0,J), 



a/3 



diag 



(4.8) 
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where v is that of s = (m ■ a) » . The exponent (matrix) v al3 is defined as 

CI g ~ s-V«\ (4.9) 

In general, V is also a matrix- valued exponent. We can easily obtain V = 1/2, J/ 2 2) , -^(33) 
for J = 0, 1/2, 1 respectively. Here is the i x i identity matrix. 

We summarize our results for the critical exponents in the table 0. There we 
show the eigenvalues of the exponent matrices v af3 and 7 a/3 , v and 7 respectively. 



Spin Factor 


Continuum Field Theory 


a 


7 


V 


V 




J=0 


free scalar field theory 


1/2 


1 





1/2 


5 (1) 


J=l/2 


Dirac field theory 


-1 


1 


1 


1 


1 + 2V2 (2) 


J=l 


massive U(l) CS field theory 


(1) 


1 


1 


1 


3 (3) 


J=3/2 













1 + V2 (4) 


J=2 




1/2 








3 (2) 


free scalar field theory 


1/2 


1 





1/2 


-3 (1) 



Table 2: Field theories and Critical exponents are obtained as scaling limit of correlation func- 
tions. The exponents 'Y and V are generally matrix-valued, and here we show only eigenvalues. 
The exponent a for J = 1 must suffer I.R. divergences. The Fisher's scaling law 7 — is(2 — r]) is sat- 
isfied. However, the hyper-scaling relation is broken for J = 1, rather inequality a > 2 — dv, (d = 3) 
is satisfied. 



5. Studies on Other Lattices 



So far we have used only the cubic lattice to define spinor fields on it. With 
another lattice structure however, one might obtain different results from those in 
the cubic lattice. Thus if we want to claim with certainty that we obtain fields of 
spin 0, \ and 1 from random walk models in three-dimensions, we have to check the 
universality of our results by taking other lattices as well. 

In this section, we mainly consider two different lattices, on which to define a 
random walk theory, a triclinic lattice and a body-centered lattice. We investigate 
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the theory on these lattices and compare the results to those on the cubic lattice. 
We show that the spin and 1/2 fields are constructed from any lattice we study 
here. On the other hand, we find an example of lattice on which the spin 1 field is 
never obtained. 

Definition First, we consider a natural definition of a random walk theory on an 
arbitrary periodic lattice L in three dimensional space R 3 . A lattice L is defined 
by putting a set of sites in R 3 periodically as a crystal with bonds which connect 
pairs of nearest neighbor sites. One must be concerned for the definition of the spin 
factor on an arbitrary lattice. 

In the case of the cubic lattice, the spin factor is determined by the oriented area 
on S 2 encircled by unit tangent vectors of a closed path. We see from the diagram in 
Fig.^| that the phase factor is defined as the area of the region on S 2 whose boundary 
is determined by points of the unit tangent vectors of a closed walk. The area of S 2 
is normalized as Air. In this way, we obtain phase factor S(/3, 7) of a path attached 
when the path turns from a direction (3 to a direction 7 by reading Fig.[| directly. 
Note that there is an ambiguity in choosing S(f3, 7) which give the same oriented 
area. 

On the basis of this spin factor, we define a random walk theory on an arbitrary 
lattice L. We assume non-backtracking walks as before. We define the amplitude of 
walks by the recursion relation 

A$ 1 (x) = ]T «Me iJ5 ^Uf 7 (x - e p ) (5.1) 

with the initial condition ( |2.3[ ). Here |eg| is the length of the bond e@ which is 
between l-th site and (I + l)-th site of the walk. S(f3,j) is the area on S 2 as 
described above. We define the free energy density with the help of A af3 as in ( |2.4| ). 

This definition of the free energy is applicable to any periodic lattice. Thus 
one might expect that spinor fields can be obtained from any lattice as well as the 
regular cubic lattice. 

In order to compute the free energy of the theory explicitly, we must require two 
additional properties on the lattice L: the length of all bonds is the same, and all 
sites of the lattice L are equivalent. 

Under these requirements, the free energy density on the lattice L can be written 
using the Fourier transformed amplitude A^{p) and the recursion relation corre- 
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sponding to Q2.5| ) : 



-(-i) 2J EEy^r(o 

1=1 a 1 



\2J 



logdet[l - kQ l (p, J)]. 



(5.2) 



Here Ql(p, J) = Ql(0, J) • -Pl(p) and a denotes the direction of a bond coming out 
of each site in the lattice. 

If we choose a lattice L, i.e., determine a relation between sites and bonds in L, 
we can fix a unit tangent vector e Q directed to a bond a (a = 1, 2, 2n). Here 2rz 
is the number of bonds come out of one site on L. We assume e a+n = —e a . Now 
we obtain the matrix -Pl(p) by an analogy with 







V 



o 

where the value p' a is defined as an inner product of p and e a , i.e. 



(5.3) 



Po 



P ■ e Q 

(Px,Py,Pz 



(5.4) 



The matrix Ql(0, J) is determined by depicting a 'diagram of the phase assignment' 
as in Fig.g. 

Note that the matrix Q(p, J) and the free energy constructed from it have the 
symmetry J — > —J as in the case of cubic lattice, i.e., Ql(p,~ J) = Ql(p,J)- 
However, the symmetry J — > J + 4 does not hold any more in general. It follows 
that the possibility of higher spin fields appearing now exists. 
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In the case of J = 0, the explicit form of Q(0, J = 0) is 

/l • • • 1 1 
. . . i . 



Q(0, J = 0) 



1 • 

1 

1 ■ 



Vi 



l l 
l l 



1 1 



1\ 



• 1 

1 

■ 1 



1/ 



(5.5) 



Some valued elements in the matrix represent the non-backtracking property of 
walks. All eigenvalues of this matrix can be shown to be either 2n — 1, 1 or —1. The 
degeneracy being 1, n and n — 1, respectively. 

The summation over all possible walks converges if and only if \k\ < k c = 
(2n — l) -1 . Since the eigenvalue 2n — 1 has no degeneracy, we expect to obtain 
the Klein-Gordon theory near this critical point as in the case of cubic lattice. The 
corresponding free energy is given by 



/(«,J = 0) 



*" d 3 p . . ri 1 - 
- — — log det 1 

-n (2tt) 3 & L 2n 



-Q L (p,J = 0)). 



(5.6) 



The singular part of the free energy (|5.6|) is 



A ((P1+P2 + ■ ■ -+P n 2 ) + 2(n - l)s) +■■■], 

(5.7) 
where 

_ 2 2n - 2 {n - l) n ~ l n n - 1 
~ (2n - l) 2 "- 1 ' 
We see the phase transition is of second order as in the case of cubic lattice. 

For other values of J, it is difficult to determine the critical behavior of the 
theory without specifying a lattice structure because the form of the matrix Q(p, J) 
depends on the lattice. 

In the following, we take two different lattices, triclinic lattice and body centered 
lattice, and examine the theory in detail especially in the case of J = 1/2 and 1. 



f L ( K ,J = 0)= TCTlog 
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Figure 3: The triclinic lattice. 




Figure 4: Definition of the phase of triclinic lattice. 

The Triclinic Lattice Here we consider the theory on a triclinic lattice, that is 
a lattice defined by inclining a cubic lattice as depicted in Fig.|3|. 

If we assume that the three bonds have the same length 1, it is characterized 
by three parameters 9, <fi and ip each of which represents the angle between two 
adjacent faces. We call such a lattice Lq^. The range of these three parameters 9, 
4> and ip is 

<9,<j),tp <7T , 7T<9 + <f) + ^<3lT. (5.8) 

Let us define the propagation matrix Qe,0,v>(P> J) f° r the lattice L 9 ^. The diagram 
of the phase assignment is given by Fig.[|, in which the area of S 2 is normalized to 
be 47T. We can read from the diagram each element of the matrix : 



20 



e 

29 



1 

—<biJ 



T -iJ 



e 2 


e (w-4>)iJ 
\ e (8-4>)iJ 



e 4>iJ 
1 

7T-2l/> 

e 2 

p (4>—ir)iJ 



iJ 







e 2 

e 2 lJ 



1 

29-7T 

e 2 

ir — 2ip 







e (7T-^)iJ 

"iJ 



-20 



e 2 



iJ 
iJ 



e (4>~ir)iJ 



e 2 lJ 

e <friJ 



-<f>ij 



-(d+4>+ir)iJ (4>+ip-2ir)iJ 



e (4>-e)iJ \ 

e {il>-<j>—ir)iJ 


e (6»+(^+7r)iJ 
e (27r-V-0)jJ 

(5.9) 



/ 



On the other hand, the matrix P(p) is defined by unit tangent vectors of three bond 
axes e' l5 e' 2 and of the lattice as 



P(P) 



V 



e w 2 











\ 



-»Pi 



(5.10) 



where = p • e' a . 

To investigate the scaling limit of the theory we must find the degeneracy of the 
eigenvalue A m of Qe,<p,i>(p, J) whose absolute value is not less than that of the other 
eigenvalues. First, we consider the J = 1/2 model. Solving the eigenvalue equation 
explicitly, we find that the six eigenvalues of the matrix Qe^^O, J = 1/2) always 
appear in degenerated pairs. We write them as Ai, Ai, A2, A2, A3 and A3. They are 
continuous functions of (9,(f>,ip) : 

Aj = \i(9,(f),ip). 

If we set 

A x (tt/2, tt/2, tt/2) > A, (tt/2, tt/2, tt/2) (j = 2, 3), 

we have the relations |Ai| > |A 2 | and |Ai| > |A 3 | for any values of 9, and ip. We 
can show this by numerical calculation. It follows that the largest absolute value of 
eigenvalues of Qe^^O, 1/2) is Ai and its degeneracy is 2. We can define the free 
energy of the theory when the diffusion constant is in the region \k\ < k c = Ai -1 . 

The free energy can be expanded in s 



fo, M (p\K, J = 1/2) = log AijPiPj + Bs 1 + 0(p\ s 3 ) 



(5.11) 



where Ay = 1,2,3) and B are real functions of (9,<f>,tp). The free energy 
fe,4>,4>(p\ K , J = 1/2) does not have rotational invariance, rather it has the symmetry 
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of an ellipsoid. This elliptic symmetry is natural in the sense that the original 
random walk on the lattice Lg^ does not have cubic but triclinic lattice symmetry 
Thus to define a rotational invariant field theory, we have to employ a new coordinate 
system which restores the rotational invariance of the theory in the scaling limit. In 
this way the free energy is rewritten in a rotational invariant form 

; w (p| K , J = 1/2) = log[V + Bs 2 + 0(p% (5.12) 

This is the same form as the cubic lattice case (c.f. (|3.3|) ). Thus, we conclude that 
we can obtain the Dirac spinor field from random walk theory also on a lattice 
by taking the scaling limit of the theory. 

Next, we consider the J = 1 case. In this case, we expect that a vector field is 
constructed on the lattice as well as on the cubic lattice. In order to obtain 

this result, it is necessary to show the triple degeneracy of eigenvalues of Q(p) whose 
absolute value is larger than that of other eigenvalues. 

Explicit calculation of the eigenvalues of Qg^^O, J = 1) however shows that 
there is no such eigenvalue for any value of 9, and ip ohter than 9 = = ip = tt/2. 

It follows that in this case we cannot obtain a vector field on the lattice Lg^ in 
general. Therefore, we conclude that the construction of the vector field from this 
random walk model is not universal. 

To see the properties of eigenvalues of Qg )t f>^(0, J) in the case of J = 1/2 and 
J = 1, we restrict the values of parameters in Qg^^(0, J) as 9 = (p = tp(= t) and 
plot the eigenvalues of Qt(0, J) (Fig.^|, 

Body Centered Lattice We consider another lattice, body centered cubic lattice. 
A part of the lattice is depicted in Fig.|^. We define a random walk on the lattice 
as a walk from the origin such that a walker standing on a site is allowed to move 
into one of the eight diagonal directions, ±e x ± e y ± e z . We require that the walker 
never pass through the bonds e x , e y and e 2 . As before, we draw the diagram of a 
phase assignment on this lattice as Fig.|8], from which the phase factor, i.e., all the 
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Eigenvalues of Q(0,J=l/2) 
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Eigenvalues of Q(0,j=l) 
6 r 




Figure 6: Eigenvalues of Q t {0, J = 1) 




Figure 7: Body centered lattice. 



24 



Figure 8: Definition of the phase of body centered lattice. 



elements of Q(0, J), of a walk can be read : 



Q(0,J) = 



V 

The matrix -P(p) is determined naturally as 

( e iP+++ 

e ip++- 

e ip+- + 



1 


gfvriJ 


1 


gfTTiJ 











1 


g-fTTiJ 


1 


g-§7riJ 


1 


g^7riJ 





1 


e" 


-nij 


1 


g|7rij 


1 




g— 7ri J 


g-§7riJ 





e~ 






1 


g — 27Ti J 


1 


g j7riJ 


g7Tl J 


g g7Tl J 










g — -| 7T'i J 


g7rij 


e -§7riJ 


1 


1 


l 


e~ 


^7ri J 


e \-KiJ 





gfTTiJ 


e ~niJ 


1 


1 


giTriJ 




1 


e -niJ 


1 





g-^TTiJ 


1 




l 


e~ 


§7™J 


1 


g^J 


giTTiJ 







1 


g§7TiJ 




1 



f(p) 







-»p++- 







(5.13) 



-ip+ — J 

(5.14) 
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Spin Factor 


Eigenvalues 


J = 


7 (1) 


1 (4) 


-1 (3) 


J = 1/2 


l + v^+v^ (2) 


l + v^-v^ (2) 


1-V3 (4) 


J = 1 


4 (3) 


(3) 


-2 (2) 


J = 3/2 


1 + (4) 


1-VQ (4) 




J = 2 


4 (2) 


2 (3) 


-2 (3) 


J = 5/2 


1-V3-V6 (2) 


1-V3 + V6 (2) 


1 + V3 (4) 


J = 3 


-5 (1) 


3 (3) 


1 (4) 



Table 3: Table of eigenvalues of the matrix Q(0, J). Numbers in () represent the degeneracy of 
the corresponding eigenvalues. 



where ... are defined as 

P±±± = ^{±Px±Py±Pz)- (5.15) 

we see from the form of the matrix Q(0, J) that the theory has the symmetry 

J -> - J, J -> J + 6. 

Using the above data, we define the free energy of the theory. The eigenvalues 
of the matrix Q(0, J) are shown in the Table |j. 

In the case of J = and 1/2, we obtain the same results as those on the cubic 
lattice at least with respect to the free energy of the theory. Namely, we obtain the 
same free energy as (|3.1|) and ( ft.3| ) up to the positive coefficient of the mass term. 

Since the matrix Q(p, J) is of order eight, fields with spin higher than 3/2 might 
come out. But the calculation of the eigenvalues of Q(p, J) show that these do not 
appear. 

We can also consider the theory on a body centered rectangular lattice. On this 
lattice, we also obtain the Dirac spinor field in the J = ~ case. Besides, the spin 1 
vector field no more appears on these deformed lattices in the J = 1 well as 

on the triclinic lattice. 
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6. Concluding Remarks 



In this paper, we construct spinor field theories in terms of a particular class of 
random walks on several three dimensional lattice structures. Since we employ such 
regularization, all physical quantities can be evaluated as finite ones at any stage 
before taking the limit of lattice spacing to be zero. We showed the convergence of 
the infinite summation over all possible walks. Therefore, the self-consistency of our 
construction is clear. 

On the regular cubic lattice, two point functions in the Klein-Gordon, the Dirac 
and the massive Chern-Simons theories are obtained for the corresponding param- 
eters J = 0, 1/2 and 1, respectively. In the J = 1/2 model there are no fermion 
doublers. In the J = 1 model three other excitations with long range correlation 
arise and govern its critical behavior, then the critical exponents do not satisfy the 
hyper-scaling relation but a corresponding inequality. This is simply due to the two 
different scale lengths. 

What is nontrivial in our results is that no spinor field theory of spin higher than 
1 is obtained. The situation is unchanged even on the body-centered lattice where 
walkers can move to eight-directions from any site. It should be studied further 
whether this ' J = 1 barrier' is a real one or an artifact of our particular choice of 
lattice. 

We have also argued the universal properties of the construction of spinor fields 
on various lattices. We see that the Klein-Gordon field is obtained on any lattice 
for J = 0. 

In the case of J = 1/2, all lattices investigated here permit the construction of 
the Dirac spinor field. A natural expectation is that the Dirac spinor field can be 
constructed on any other lattice. 

The authors thank T. Hara, K. R. Ito, N. Sakai, G. Semenoff and H. Tasaki for 
important comments and discussions. One of them (S.K.) thanks menbers of Yukawa 
Institute for Theoretical Physics for kind hospitality extended to him. They are 
grateful to J. Ambj0rn, P. Kurzepa and P. Orland for notifying them of references. 
They thank P. Crehan for reading the manuscript. This work is supported in part by 
Grant-in- Aid for Scientific Research (S.K.) from the Ministry of Education, Science 
and Culture. 
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Appendix 



We prove the convergence property of the free energy defined by the first equation 
of the eq. Q2.7p . We show that the convergence radius of the infinite sum as a power 
series in k 



OO 1 

A(k,J) = txJ2j^Q(p,J)} 1 
1=1 1 



(A.l) 



is \k c \, where n c 1 is the eigenvalue of Q(0, J) with the largest absolute value. Let 
U(p, J) be a unitary matrix which triangulates the matrix Q(p, J) as 



/Ai(p,J) 



u( P ,j)Q(p,j)u(p, jy 



A a (p, J) 



V 







(A.2) 



A 6 (P, J)/ 



Note that the diagonal elements Aj(p, J) (i = 1, 6) are identical to the eigenvalues 
of the matrix Q(p, J). By using this expression, we can deform the sum ( |A.1| ) as 
follows : 



A(k, J) 



tr jU(p, jy [kU(p, J)Q(p, J)U(p, J)t U(p, J) 
i=i 1 



OO 1 

1=1 1 



/Ai( P , J) 



K 



\ 



V 







Up, J) J 



6 L 2 

lim J^yf/cA^p, J)]'. 

L ^°° i=i ;=i « 



(A.3) 



If |«;Aj(p, J) | < 1 (i — 1, 6) holds, this summation converges to be 



^log(l - «Ai(p, J)), 



i=i 



(A.4) 
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which is also represented as 



logdet [1 — kQ(p, J)] . (A. 5) 

Thus the free energy converges if k is in the region \k\ < \k c \ where n c = k c {J) 
is the inverse of the eigenvalue whose absolute value is the largest among the set of 
eigenvalues of matrices Q(p, J) for a fixed J. 

Next, we show that the eigenvalue with the largest absolute value for a fixed J 
is given by that for p = 0. Let Aj {i = 1, 6) be the eigenvalues of Q(J, 0). We set 
that the absolute value of Ai (or Xq) is less (or more) than or equal to any other Aj. 
Given a unit vector x, there exists a unit vector y = (y 1: • ■ ■ , y 6 ) such that 

xt Q( J, 0)tQ(J, 0) x = ]T I^IV- (A.6) 

i=i 

It follows that 

Ai 2 <x f g(J, 0YQ{J, 0)x< A 6 2 . (A. 7) 

For any eigenvalue a of Q( J, p) = Q(J, O)P(p), we can choose a unit eigenvector 
x : Q(J, 0)P(p)x = ax . Then, 

(Px)tQ 2 Px = \a\ 2 . (A.8) 

By applying the formula ( |A.7|) to ( |A.8| ), we obtain 

Ai 2 < \a\ 2 < A 6 2 . (A.9) 

It has been shown that the eigenvalue of the largest absolute value is that of 
p = 0. Therefore the infinite sum eq.( |A.l| ) converges uniformly with respect to p 
for \k\ < \k c \. 
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